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ABSTRACT 



■ Received jdate/,; Accepted jdate^ 

o 

I ■ Aims. We investigate the environment of the infrared dust bubble N65 and search for 

2 ■ evidence of triggered star formation in its surroundings. 

C/5 , Methods. We performed a multiwavelength study of the region around N65 with data 

^ ' taken from large-scale surveys: Two Micron AU Sky Survey, GLIMPSE, MIPSGAL, 

SCUBA, and GRS. We analyzed the distribution of the molecular gas and dust in the 
^ ' environment of N65 and performed infrared photometry and spectral analysis of point 

sources to search for young stellar objects and identify the ionizing star candidates. 



in 



OO ' Results. We found a molecular cloud that appears to be fragmented into smaller clumps 

along the N65 PDR. This indicates that the so-called collect and collapse process may 



O 



be occurring. Several young stellar objects are distributed among the molecular clumps. 
, They may represent a second generation of stars whose formation was triggered by the 

bubble expanding into the molecular gas. We identified 0-type stars inside N65, which 
are the most reliable ionizing star candidates. 



X 
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1. Introduction 

Since the early work of Elmegreen & Lada (1977), it is known that during the supersonic 
expansion of an HII region, a dense layer of material can be collected between the ionization 
and the shock fronts. This layer can be fragmented into massive condensations that may 
then collapse to form new massive stars and/or clusters. Several observational studies have 
inferred that this process, known as "collect and collapse" , triggers massive star formation 
(see e.g., Pomares et al. 2009; Zavagno et al. 2007, and references therein). As Deharveng 
et al. (2005) point out, only the presence of either a dense molecular shell surrounding the 
ionized gas of an HII region or massive fragments regularly spaced along the ionization 
front, can prove that we arc dealing with the collect and collapse process. 

Using Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) data, 
Churchwell et al. (2006, 2007) cataloged almost 600 infrared (IR) dust bubbles, full or 
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partial rings bordered by a photodissociation region (PDR) and detected mainly at 8 /im, 
which usually enclose ionized gas and hot dust observed at 24 ^m. Most of these bubbles are 
HII regions and several have morphologies that are indicative of triggered star formation. 

On the other hand, Cyganowski et al. (2008) also using GLIMPSE data, identified more 
than 300 extended 4.5 fim sources, so-called "extended green objects (EGOs)", applying 
the common coding of the [4.5] band as green in three-color composite Infrared Array 
Camera images. According to the authors, an EGO is probably a massive young stellar 
object (MYSO) driving outflows. The extended emission in the 4.5 /im band is produced 
by H2 {ly = 0-0, 8(9,10,11)) lines and CO (i^ = 1 — 0) band heads, which are excited by 
the shock of the outflows propagating in the interstellar medium (ISM). The majority of 
EGOs are associated with infrared dark clouds (IRDCs) and some lie over the border of 
IR dust bubbles. 

From the catalog of Churchwell et al. (2006), we select the IR dust bubble N65, which 
harbors the EGO G35. 03-1-0. 35 in one of its borders (indicative of active star formation), 
to study its environment. 

N65 is a complete (i.e., closed ring) IR dust bubble centered on a2ooo — 18''54™02*'. 
(52000 = +01°59'30" {I = 35?000, b = +0?332) with a radius of about 2'.6 (Churchwell 
et al. 2006). Close to N65, over one of its borders, centered on 0:2000 = 18^54'°04.2'^, 
(^2000 = +02°01'33.9" {1 = 35?035, b ^ -h0?338) lies the IR source IRAS 18515+0157, 
where several molecular lines tracers of star formation were detected. Maser emission from 
II2O and OH was detected by Forster & Caswell (1989), as well as methanol (Caswell 
et al. 1995). Bronfman et al. (1996) and Jijina et al. (1999) also detected CS and NH3, 
respectively. The detection of formaldehyde at v^sR ^ 57 km s~^ was used by Watson 
et al. (2003) to determine that this star-forming region lies at the kinematic distance of 10 
kpc, the farthest possible given the distance ambiguity of the first Galactic quadrant. On 
the other hand, close to this region, at I = 35?015, b = -|-0?356 and also across the border 
of N65, Anderson et al. (2009) cataloged the UCHII region G35.02-K0.35 with a velocity of 
57.2 km s~^, which they locate at the near distance of ~ 3.6 kpc. 

In this work, we present a molecular and IR study of the environment surrounding 
the IR dust bubble N65 to explore its surrounding ISM and search for signatures of star 
formation. We describe the data used in Sect. 2, the results and the discussion of them are 
presented in Sect. 3 and Sect. 4 summarizes the results. 

2. Data 

We analyzed data extracted from five large-scale surveys: the Two Micron All Sky 
Survey (2MASS)\ Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE), 
MIPSGAL, GRS^, and the SCUBA Legacy Catalogue. GLIMPSE is a mid-infrared survey 
of the inner Galaxy performed with the Spitzer Space Telescope. We used the mosaicked 

^ 2MASS is a joint project of the University of Massachusetts and the Infrared Processing and 

Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space 

Administration and the National Science Foundation. 
^ Galactic Ring Survey (Jackson et al. 2006) 
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images from GLIMPSE and the GLIMPSE Point-Source Catalog (GPSC) in the Spitzer- 
IRAC (3.6, 4.5, 5.8, and 8 fim). IRAC has an angular resolution of between I'.'S and 1'.'9 
(see Fazio et al. 2004 and Werner et al. 2004). MIPSGAL is a survey of the same region as 
GLIMPSE, using the MIPS instrument (24 and 70 ^.m) on Spitzer. The MIPSGAL resolu- 
tion is 6" at 24 jum. The GRS is being performed by the Boston University and the Five 
College Radio Astronomy Observatory (FCRAO). The survey maps the galactic Ring in 
the ^^CO J=l-0 line with an angular and spectral resolution of 46" and 0.2 km s^-'^, respec- 
tively (see Jackson et al. 2006). The observations were performed in both position-switching 
and on-the-fly mapping modes, achieving an angular sampling of 22". The SCUBA Legacy 
Catalogues'^ are two sets of continuum maps and catalogs using data at 450 and 850 fim 
obtained with the Submillimetre Common User Bolometer Array (SCUBA) with angular 
resolutions of 19" and 11", respectively (Di Francesco et al. 2008). 



3. Results and discussion 

Figure 1 (left) shows a Spitzer-IRAC three color image of N65. The three IR bands are 3.6 
/im (in blue), 4.5 ^m (in green), and 8 fim (in red). The location of the EGO G35.03+0.35 
(Cyganowski et al. 2008) is indicated. Figure 1 (right) displays a composite two-color image 
towards N65. Red and green represent the Spitzer-IRAC emission at 8 /im and the Spitzer- 
MIPSGAL emission at 24 /im, respectively. Both figures clearly show the PDR visible in the 
8 /im emission, which originates mainly in the polycyclic aromatic hydrocarbons (PAHs). 
The PAH emission delineates the HII region boundaries because these large molecules are 
destroyed inside the ionized region, but are excited in the PDR by the radiation leaking 
from the HII region (Pomares et al. 2009). The 24 /im emission, displayed in green in Fig. 
1 (right), corresponds to hot dust, distributed mainly towards the eastern border of N65 
that coincides with the central position of IRAS 18515-1-0157. 




35,120 35 080 36 040 35 000 34 960 35 120 35 030 35 040 35OO0 34 960 

Galactic longitutie Galaclic longitutJe 

Fig. 1. Mid-IR emission of the IR dust bubble N65. Left: Spitzer-IKAC three-color image (3.6 /im 
= blue, 4.5 /im = green, and 8 fim = red). The location of the EGO G35. 03-1-0. 35 is indicated. 
Right: red is the Spitzer-IRAC 8 /im emission and green is the Spitzer-MIPSGAL emission at 24 
/im. 

^ http:/ /www2.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/scubalegacy/ 
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From Fig. 1, a probably smaller IR dust bubble can be discerned in the southwest border 
of N65, centered on I = 34?96, b ~ +0?31, which is displayed in Fig. 2. The emission at 
8 ^m (Fig. 2 left) shows an almost circular FDR, which encloses a region of about f 8" in 
diameter. At the center of the bubble, a source is visible with emission in the 3.5, 4.5, 8, 
and 24 fim bands, which is probably the exciting star of this structure surrounded by hot 
dust emitting at 24 fim (green in Fig. 2 right). 




Fig. 2. Mid-IR emission of the possible smaller IR dust bubble that lies in the southwest border 
of N65. Left: Spitzer-IRAC three-color image (3.6 fj,m = blue, 4.5 /j,m — green and 8 fj,m — red). 
Right: red is the Spitzer-IRAC 8 /im emission and green is the Spitzer-MIPSGAL emission at 24 
/im. 



3.1. Molecular and dust analysis 

Wc inspected the molecular gas around N65 from the GRS data in the whole velocity range 
and foimd an interesting feature around v ^ 50 km s^^. Figure 3 shows the integrated 
velocity maps of the ^'^CO J=f-0 emission every 1.05 km s^^ between 47 and 56 km s^^. 
The circle shows the position and size of N65. A molecular cloud with an arc-like shape 
open to the south is evident that encloses the IR bubble over one of its borders. The 
molecular gas exhibits several clumps along the PDR. The distribution and morphology of 
this material suggests that the collect and collapse process might be occurring. 

In Fig. 4, we compare the ^'^CO emission integrated within the interval 47 and 55 km 
s~^ with the continuum dust emission at 850 /xm extracted from the SCUBA (in yellow 
contours). In this figure, it is evident that the dust emission coincides with the brightest 
molecular clump centered on I ^ 35?02, b ^ 0?35 and with the position of the EGO 
G35. 03-1-0. 35. The presence of the EGO together with all the previous results mentioned 
in Sect. 1 provide clear evidence that this region might be active in star formation. 

If we assume that the molecular gas in the velocity range considered before is physically 
associated with N65, we can adopt ^ 50 km s^^ as the systemic velocity of the infrared dust 
bubble. According to the galactic rotation model of Fich ct al. (1989) (with Rq — 8.5 kpc 
and Vq — 220 km s~^), we obtain kinematic distances of cither 3.5 or 10.5 kpc. Considering 



A. Petriella et al.: The 



environment of the infrared dust bubble N65 



5 



2 4 6 




35.16 35.08 35.00 34.92 



Galactic Long. 

Fig. 3. Integrated velocity maps of the ^'^CO J=l-0 emission every 1.05 km s~^. The grey scale is 
in K km s~^. The contour levels are 2, 4 and 6 K km s~^. The circle represents the position and 
size of N65. 



that the detection rate for IR dust bubbles in GLIMPSE peaks at the distance of 4.2 kpc 
within a horizon of 8 kpc (Churchwell et al. 2006), the bubbles at the near kinematic 
distance are more likely to be detected than those at the far distance. Thus, we assume a 
distance of 3.5 kpc for N65, in coincidence with the UCHII region G35.02+0.35 cataloged 
by Anderson et al. (2009). 

Using the ^'^CO J=l-0 line and assuming local thermodynamic equilibrium (LTE), we 
estimate the H2 column density toward the brightest molecular clump shown in Fig. 4. We 
use 

NC^CO) ^ 2.42 X 1014 ^'iZt/T V 
1 - exp(~5.29/Tc^) 

to obtain the ^'^CO column density, where T13 is the optical depth of the line and following 
Anderson et al. (2009) we considered T^x = 20 K. We assume that the ^'^CO emission is 
optically thin and use the relation N(H2)/N(i3co)-- 5 x 10^ (e.g., Simon et al. 2001) to 
estimate a column density of N(H2) ^ 1.6 x 10^^ cm^^. From this value, we estimate the 
mass of the molecular clump to be ~ 2 x 10"^ M0. This value was obtained from 

M = ^mH;^[i5' f^N(H2)], 

where is the solid angle subtended by the ^'^CO J=:l-0 beam size, is the hydrogen 
mass, /i is the mean molecular weight assumed to be 2.8 by taking into account a relative 
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Fig. 4. "CO emission integrated between 47 and 55 km s , the contours levels are 19, 23, and 
33 K km s^^. The yellow contours correspond to the 850 fim continuum emission obtained from 
SCUBA and the levels are 0.3, 0.6, 0.9, and 2.3 Jy beam~\ The beams of each emission are 
included in the bottom left corner. 

helium abundance of 25 %, and D is the distance. Summation was performed over all the 
observed positions within the 33 K km s~^ contour level (see Fig. 4). Finally from this 
mass value, we derive a density of n ^ 10^ cm"'^. 



Using the 850 /im continuum emission, we estimate the dust mass from the relation of 
Tej et al. (2006): 



where S'^ is the flux density at the frequency v. We use S353GHZ — 19.42 Jy as obtained 
by Di Francesco et al. (2008) for this region. We assume that T^i the dust temperature is 
20 K, /? the dust emissivity index is 2.6 for the assumed dust temperature in this region 
according to Hill et al. (2006), and D the distance to N65 is 3.5 kpc. This equation assumes 
a standard opacity of ki2ooGHz =0.1 cm^ g^^. We obtain a dust mass of Mdust ~ 60 M©. 
The dust-to-gas mass ratio is then ~ 0.03, which is 3 times higher than the canonical 
Galactic ratio. 

3.2. IR photometry 

We analyze the distribution of IR point sources in the surroundings of N65 to search for 
signs of star formation and try to identify the ionizing-star candidates. Figure 5 shows 
the [5.8]-[8.0] versus [3.6]-[4.5] color-color (CC) diagram for the sources extracted from 
the GLIMPSE Point Source Catalog in the Spitzer-IRAC bands within a circle of 8' in 
radius centered on I = 35? 00, b = +0?32. The size of this region completely covers the 
extension of N65 and the surrounding molecular clumps, where stars may be forming. We 
only consider sources with detection in the four bands. The regions in the figure indicate 
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the stellar evolutionary stages based on the criteria of Allen et al. (2004) : class I sources are 
protostars with circumstellar envelopes, class II are disk-dominated objects, and class III 
are main sequence and giant stars. We search for ionizing star candidates among class III 
sources, while class I and II sources are chosen to be the YSO candidates. In the following 
sections, we report the results of these analyses. 
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Fig. 5. GLIMPSE color-color diagram [5.8]-[8.0] versus [3.6]-[4.5] for sources within a circle of 
8' in radius centered at N65. We only consider sources with detection in the four Spitzer-YRAC 
bands. The regions indicate the stellar evolutionary stage as defined by Allen et al. (2004). Class 
I sources are protostar with circumstellar envelopes; class II are disk dominated objects; and class 
III are main sequence and giant stars. 



3.3. Identifying tlie exciting star(s) of N65 

We attempt identify or at least propose the possible exciting stars(s) of N65. 

According to Churchwcll ct al. (2006), IR dust bubbles may be produced by O- and/or 
B-type stars. Given that we do not find any cataloged massive star toward N65, we look 
for the exciting star(s) among the class III sources because they include main sequence 
stars. The ionizing star(s) arc not always located close to the geometrical center of the 
bubble because of their proper motion. However, we expect to find them inside the PDR 
as IR dust bubbles are young objects and the ionizing star(s) may not have enough time 
to overcome the shell. To perform near IR photometry using the 2MASS JHK bands, we 
select class III objects within the N65 PDR that have detections in the aforementioned 
near IR bands. Figure 6 shows the {H-K) versus {J-H) color-color (CC) diagram of the 
22 sources found using this criteria. Figure 7 shows the K versus {H-K) color-magnitude 
(CM) diagram. Finally, in Fig. 8 we show the location of these sources. 

To identify possible exciting star(s) candidates among the 22 sources found, we select 
objects that meet both of the following criteria simultaneously: (i) according to their posi- 
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tion in the CC diagram (Fig. 6), the sources may be main sequence stars and (ii) according 
to their position in the CM diagram (Fig. 7), the sources may be O- or B-typc stars. These 
sources are indicated by red numbers in Figs. 6 and 7 and red crosses in Fig. 8. 



y» M .22 



:.'",^a'„io 



Fig. 6. (H-K) versus (J-H) color-color (CC) diagram of the class III objects within the N65 PDR. 
The two solid curves represent the location of the main sequence (thin line) and the giant stars 
(thicker line) derived from Bessell & Brett (1988). The parallel dashed lines are reddening vectors. 
We have assumed the interstellar reddening law of Rieke & Lebofsky (1985) {A j / Ay =0.282; 
Ah/Av=0.175 and Ak/Av=0.112). The sources with red numbers are the ones selected to be the 
candidates for exciting stars of N65. 

As an additional constraint of the ionizing star(s) candidate, we derive the spectral 
energy distribution (SED) of the selected sources by fitting the fluxes in the IRAC and 
2MASS bands with a Kurucz photospheric model using the tool developed by Robitaille 
et al. (2007) and available online'*. The fitting tool requires a distance range and a range 
of visual extinction. We assume a distance of 3.5 kpc and derive the extinction for each 
source from the 2MASS CC diagram (H-K) versus (J-H) presented in Fig. 6. Once the 
SED is obtained, we check whether the fitted effective temperature is consistent with the 
temperature of an O- and B-typc star. Our analysis is only qualitative because of the 
uncertainty in the determination of the spectral type from the CM diagram and because 
the fitting photosphere model extends only to 50,000 K. 

The fitted effective temperatures of stars 5, 10, and 16 obtained from the SED is ^ 
6,000 K and for star 20 is 20,000 K. These values are too low for a O- or early B-type 
star. For the remaining sources 8, 9, 12, 15, 17, 18, 20, and 21 the fitted temperature is ^ 
50,000 K, which is in closer agreement with the effective temperature expected for an O- 
type star (see Schaerer & de Koter 1997). Thus, based on their position inside the bubble, 
we suggest that sources 8, 12, and 15 are the most reliable ionizing-star candidates of N65. 



* http;/ /caravan. astro.wisc.edu/protostars/ 
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Fig. 7. (H) versus (H-K) color-magnitude (CM) diagram of the class III objects within the N65 
PDR. The position of the zero-age-main-sequence stars with A,j = mag at a distance of 3.5 kpc 
are indicated by black stars. The reddening curve for an 03 star is shown with a black line. The 
sources with red numbers represent are the ones selected to be candidate exciting stars of N65. 
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Fig. 8. Possible exciting star(s) over the 8 /xm emission of N65. The red crosses are the sources 
selected, according to the CM and the CC diagrams, to be candidate exciting stars. 



3.4. Star formation around N65 

Figure 9 shows the distribution of both class I (yellow crosses) and class II (magenta 
crosses) point sources around N65. We note the presence of a group of these sources right 
upon the molecular clumps that surround the bubble above its PDR. As molecular clouds 
are star birth places, we expect some of these sources to be YSOs whose formation could 
have been triggered by the expanding IR dust bubble. 

We perform a fitting of the fluxes of YSO candidates in the IRAC and 2MASS bands, to 
derive their spectral energy distribution (SED) and constrain their evolutionary stage. We 
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Fig. 9. Two-color imag e towards N65: "CO emission integrated between 47 and 55 km s ^(in 
green with white contours) and 8 /im emission (in red) . We show the distribution of class I (yellow 
crosses) and class II (magenta crosses) sources around N65 derived from the CC diagram of Fig. 
5. We have labeled the sources located upon the molecular cloud and close to the PDR. 



limit our study to the sources superimposed on the molecular cloud around N65 (see Fig. 
9). The SED wa.s obtained using the tool developed by Robitaille et al. (2007) available 
online^. This tool consists of a grid of precomputed radiative transfer models with a fast 
minimization algorithm. We select the models that satisfied the condition 

- xLn < '^N, (1) 

where Xmin minimum value of the among all models, and A'^ is the number of 

input data fluxes (fluxes specified as upper limit do not contribute to N). Hereafter, we 
refer to models satisfying Eq. (1) as "selected models". 

To relate the SED to the evolutionary stage of the YSO, Robitaille et al. (2006) de- 
fined three different stages as a complementary classification of Allen's classes. The stage 
classification is based on the values of the central source mass M^,, the disk mass Mdisk, 
the envelope mass M^nv, and the envelope accretion rate Menu of the YSO. Stage I YSOs 
are those that have Menv/AIi, > 10^^ y^^, i-c, protostars with large accretion envelopes; 
stage II are those with Mdisk/ Mi, > 10^^ yr^^ and Me„i,/A/* < 10^^ yr^^ , i-e., young 
objects with prominent disks; and stage III are those with Mdisk < 10^^ yr"-^ and 
Menv/Mi, < 10~^ yr~^, i.e., evolved sources where the flux is dominated by the central 
source. The stage classification is based on the physical properties of the YSO, rather than 
^ http;/ /caravan. astro.wisc.edu/protostars/ 
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on the spectral index derived from the slope of its SED. However, as pointed by Robitaille 
et al. (2006), the stage classification is equivalent to the class classification for a great num- 
ber of models. Therefore, an easy way of identifying young objects is looking for sources 
that can be fitted only by selected models with Menv > 0. Sources that are fitted by mod- 
els with zero and non-zero values of M^^v may be more evolved YSOs whose evolutionary 
stage is not constrained by the data (see Poulton et al. 2008). 

In Table 1, we report the main results of the fitting output for class I and II sources 
projected onto the molecular cloud. The YSO identification (Col. 1) is taken from Fig. 9. In 
Col. 2 and 3, we report the pcr data point of the best- fit model and the number on models 
satisfying Eq. (1), respectively. The remaining columns report the physical parameters of 
the source, specifying the range of values of the selected models: central source mass, disk 
mass, envelope mass, and envelope accretion rate, respectively. These results permit us to 
confirm the presence of young objects in the vicinity of the IR dust bubble N65. 



Table 1. Parameters of class I and II YSO candidates derived from the SED fitting. 



Source 


'X.min / ^ 


n 
















(Mo) 


(Me) 


(Mo) 


(Mo/yr) 


1 


0.13 


484 


1 - 11 


8.3 X 10"** - 0.4 


6.4 X 10"^ - 1800 


- 0.003 


2 


2.2 


41 


5-8 


7.1 X 10""* -0.7 


1.1 X 10"** - 1200 


- 0.002 


3 


0.1 


84 


1-11 


1.1 X 10"*^ - 0.71 


1.7 X 10"* - 140 


- 0.001 


4 


0.2 


273 


0.2 - 11 


5.5 X 10"''^ -0.6 


3.8 X 10"^ - 360 


- 0.003 


5 


5 


4 


9 





28 


3.8 X 10"^ 


6 


2.2 


2 


4.6 


0.014 


37 


2.1 X 10~^ 


7 


10 


9 


4-5 


7.4 X 10"^ - 0.02 


0.1 - 4 


5 X 10"^ - 9 X lO"'^ 


8 


1.3 


1 


9.6 


0.75 


130 


2.1 X 10"" 


9 


5 


61 


2-7 


1.6 X 10"'' - 0.13 


1.8 X 10"^ - 17 


- 5.4 X 10"* 


10 


4 


9 


5 


4.5 X 10"^ - 0.001 


0.09 - 1.5 


9.5 X 10"^ - 1.4 X lO"'' 


11 


3.3 


10 


2-4 


6.5 X 10"^ -0.03 


0.02 - 16 


4 X 10"'^ - 5 X 10"* 


12 


11 


2 


3-4 


0.016-0.02 


2-8 


3.5 X 10"® - 2 X 10"^ 


13 


2.4 


13 


1 - 7 


1.1 X 10"*^ - 0.3 


0.006 - 43 


3.5 X 10"* - 1.1 X 10"* 


14 


1.4 


298 


1-11 


3 X 10"'' -0.7 


5 X 10"^ - 140 


- 0.002 


15 


1 


3 


8-9 


0.03-0.2 


70 - 140 


1.4 X 10"* - 2 X 10"* 


16 


0.03 


430 


1 - 8 


3 X 10"* - 0.4 


1 X 10"® - 23 


- 4 X 10"* 


17 


3 


2 


4.6 


0.014 


37 


2.1 X 10"* 


18 


0.8 


13 


2-9 


6.7 X 10"'' -0.7 


0.1 - 140 


3 X 10"*^ - 4 X 10"* 


19 


0.4 


85 


2-21 


0- 1.1 


0.4 - 1300 


4.1 X lO"'"^ - 0.002 


20 


2 


2 


9-10 


0.3 - 0.8 


65 - 130 


2 X 10~* 


21 


0.003 


1652 


0.7-15 


6 X 10"'^ -0.9 


4.5 X 10"^ - 560 


0-0.001 


22 


1.3 


50 


3-15 


9 X 10"'-* - 0.2 


1.4 X 10"* - 200 


- 6 X 10"* 



Wc identify 12 sources that are only fitted by models with Menv > 0; they are sources 
5, 6, 7, 10, 11, 12, 13, 15, 17, 18, 19, and 20. These are presumably young objects embedded 
in prominent envelopes. In every case, the selected models are prevalently stage I with a 
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few stage II. The only exception is source 13, which is fitted by both stage I and stage III 
models. 

Sources 5, 15, and 20 appear to be the youngest YSOs. They may all be massive central 
sources (between 8 and 10 Mq) surrounded by envelopes of several solar masses. In the 
case of 15 and 20, their ages are estimated to be ^ 10"* years. Regarding the source 5, its 
position is coincident with the ultra-compact HII region G35.02+0.35 and its age is about 
10'^ years. This is a really young massive source and a clear sign that star formation is still 
occurring in the vicinity of N65. 

The selected models (i.e., models satisfying Eq. 1) for sources 7, 10, 11, and 12 spread 
over a similar range of values among them. These YSOs contain a central source of mass 
less than 5 Mq and ages between 10'* and 10^ years. 

For source 19, we obtain a large number of selected models (85) because for this source 
we only fitted the four IRAC fluxes. However, all of them are stage I so this source is 
probably a young object. The large spread in the range of parameters does not allow us to 
constraint its mass and age. 

Projected onto the molecular cloud, we also identify sources with zero and non-zero 
values for Menv- sources 1, 2, 3, 4, 9, 14, 16, 21, and 22. The evolutionary stage of these 
sources cannot be assured without doubt because the parameters of the selected models 
are spread over a wide range. We should include fluxes measured at other wavelengths to 
reduce the number of models satisfying the condition deflned in Eq. (1). However, we can 
identify sources that are in their early stages of evolution (embedded or disk dominated) 
because the selected models are either stage I or II. This is the case for sources 1, 2, 4, 14, 
and 16. For the remaining sources, 3, 9, 16, and 22, the selected models are distributed 
between stages I, II, and HI so their evolutionary stage cannot be established without 
controversy. 

We note that source 21 is located at the center of the smaller IR bubble introduced in 
Sect. 3 that lies right upon the PDR of N65 (see Fig. 2). We fitted the SED of the source 
and obtained more than 1600 models satisfying Eq. (1) and, consequently, a wide range of 
physical parameters for the YSO, so we cannot establish its evolutionary stage. However, 
its location in the center of the smaller IR dust bubbles suggests that this source may be a 
second generation massive star blowing a stellar wind strong enough to evacuate a bubble 
around itself. More observations are needed to investigate this scenario. 

In conclusion, the SED study confirms the existence of several YSOs possibly embedded 
in the molecular clumps that surround N65, which is indicative of the collect and collapse 
process occurring in this region. 

3.5. SED analysis of the EGO G35.03+0.35 

Wc derive the SED for the EGO G35.05+0.35 using the fiuxes in the Spitzer-IRAG 4.5, 
5.8, and 8.0 /im bands (the GPSC reports no detection in the 3.6 yum band and the source 
is absent in the 2MASS catalog) . We also use fiuxes in the far IR of the SCUBA bands of 
450 and 850 /im and in the band of 1200 ^m measured by the SIMBA bolometer at SEST 
(Hill et al. 2006). Both the SCUBA and SIMBA datasets have a lower angular resolution 
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than that of the GPSC and the measured fluxes may include contributions from other 
sources near the EGO. For this reason, we perform a fitting of the SED by considering 
two different cases: a) assuming that the far IR data represent the flux emitted only by 
the EGO, the SED of the best-fit model being shown in Fig. 10 (left); and b) using the 
far IR data as an upper limit (i.e., assuming that the emission originates in the EGO and, 
probably, other sources around it), the SED of the best- fit model in this case being shown 
in Fig. 10 (right). In Table 2, we report the physical parameters of the selected models, 
considering both SCUBA and SIMBA data as flux emitted only by the EGO (case a) and 
upper limits (case b). The column designation is the same as in Table 1. 




X (|im) X (|im) 



Fig. 10. Best-fit SED model for the EGO G35.05-f 0.35. The points in blue are the measured 
fluxes at the Spitzer-IKAC bands of 4.5, 5.8, and 8.0 /im, the SCUBA bands of 450 and 850 /^m, 
and the band of 1,200 nm measured with the SIMBA bolometer at SEST. Left: SCUBA and 
SIMBA data are assumed to represent flux emitted only by the EGO. The total flux (black line) 
is completely dominated by the flux from the envelope as would be expected in a young stellar 
object at its earlier stages of evolution. Right: SCUBA and SIMBA data are taken to be an upper 
limit. The total flux at the longer wavelengths is dominated by the envelope, while at the shorter 
wavelengths the main contributions originates in the disk (red line) and the central source (blue 
line) . 



Table 2. Physical parameters of EGO G35. 05-1-0. 35 derived from the SED fitting for the selected 
models satisfying Eq. (1). 



Case 


Xmin / ^ 


n 


M* 














(Mo) 


(Mo) 


(Mo) 


(Mo/yr) 


a 


30 


2 


27-33 


0-0.3 


1900 - 2500 


0.004 - 0.006 


b 


3.3 


16 


6-10 


0.005 - 0.8 


0.7-1100 


2 X 10"^ - 0.002 



According to Fig. 10 [left), if we consider fiuxes at the longer wavelengths to be emitted 
only by the EGO, the total fiux of the best-fit model is completely dominated by the 
envelope. As pointed out by Robitaille et al. (2006), the envelope of sources with high 
accretion rates becomes optically thick. In our case of Me„„ '~ lO^'^ Mq /yr, radiation 
emitted by both the central source and the disk is probably absorbed and re-emitted by 
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the envelope of gas and dust. As can be seen in Table 2, the physical parameters of the EGO 
indicate that it is a massive central source surrounded by a massive envelope. According 
to Fig. 10 (right), if we take SCUBA and SIMBA measurements to be upper limits, the 
total flux of the best-fit model is still dominated by the envelope but only at the longer 
wavelengths. In the near and mid IR, the main contribution to the total flux comes from 
the disk. The accretion rate values that we derived for the envelope are lower than those 
derived in the previous case (see Table 2). Thus, the envelope is expected to be more 
transparent to radiation coming from the inner regions (disk and central source). 

Both cases show that the fitted SED of the EGO is indicative of a massive central source 
surrounded by a massive envelope. The age of the YSO is estimated in between 10"* and 
10^ years. From these results, we suggest that the EGO G35.05-I-0.35 is a candidate massive 
young stellar object at the earlier stages of evolution with outflowing activity, in agreement 
with that proposed by Cyganowski et al. (2008). The presence of this object is a clear sign 
that star formation is occurring in the vicinity of N65. 

4. Summary 

Wc have conducted a study of the surroundings of the IR dust bubble N65 using large- 
scale surveys and archival data in the radio and infrared bands. The main results can be 
summarized as follows: 

1. We discovered the presence of a molecular cloud with an arc-like shape enclosing N65 
over its east, north and west borders and centered on a velocity of 50 km s^^. The 
cloud appears to be fragmented into smaller molecular clumps that arc seen right upon 
the PDR. Based on the distribution and morphology of the molecular gas between 47 
and 55 km s^^, we suggest that the collect and collapse process may be occurring. 

2. The EGO G35.05-f0.35 lies in the most intense molecular clump, which belongs to the 
discovered molecular shell. From the ^■^CO J=l-0 emission, wc estimated its molecular 
mass and density to be ~ 2 x 10"^ Mq and ^ 10^ cm^^^, respectively. We found that 
the dust contimmm emission, traced by the 850 fiia band, also peaks in this region. We 
estimated a dust-to-gas ratio of ~ 0.003 for this clump. 

3. From the analysis of the infrared photometry of point sources around N65, we identified 
three O- or early B-type stars as the most probable candidates responsible for the 
creation of N65. In addition, we found many YSOs distributed right upon the molecular 
clumps around N65. The analysis of the SED of these sources confirms that many of 
them are in the embedded or disk-dominated stage. This suggests that the gas collected 
into the molecular clumps may have collapsed to form a second generation of stars, 
showing another galactic case of the collect and collapse process. 

4. The SED of EGO G35. 05+0.35 as derived by fitting the fluxes at the mid and far IR, 
confirms that it is a massive young stellar object. 
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